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ABSTRACT 

The  thermal  desorption  of  Xe  from  the  (110)  plane  of  W  was  investigated. 
Desorption  could  be  followed  layer  by  layer  for  3  adsorbed  layers,  each  con¬ 
taining  equal  amounts  of  Xe.  In  each  case  desorption  rates  were  nearly  coverage 
independent  (i.e.  of  zero  order)  until  0  =  0.25*0.3,  (in  the  layer  under  study) 
at  which  point  a  sharp  transition  to  first  order  desorption  occurred.  If  only 
1  monolayer  of  Xe  was  adsorbed  small  amounts  of  preadsorbed  oxygen  led  to  first 
order  desorption  over  the  entire  coverage  regime.  Possible  mechanisms  for  the 
zero  order  regime  in  terms  of  desorption  from  the  edges  of  islands  with  con¬ 
stant  total  island  perimeter,  or  via  a  more  complicated  mechanism  are  given. 

The  activation  energies  of  desorption  are  shown  to  be  less  than  the  enthalpies 
of  desorption  for  all  layers  and  the  meaning  of  this  result  is  discussed.  Re¬ 
sults  for  Xe  adsorbed  on  0/W  layers  preheated  to  9QK  and  900K  are  also  presented 
INTRODUCTION 

The  desorption  of  inert  gases  from  metal  surfaces  would  seem  to  present  a 
particularly  simple  example  of  desorption  kinetics  since  the  complications  of 
dissociation  on  adsorption  and  recombination  on  desorption  are  absent.  One 
example  in  the  literature,  the  study  of  Erikson  and  Yates[l]  on  the  desorption 
of  Xe  from  the  W(lll)  surface  in  fact  seems  to  bear  this  out.  Desorption  was 
found  to  be  first  order  with  an  activation  energy  of  9*3  Kcal/mole  and  a  fre- 
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quency  factor  of  ~10^sec  Some  indication  of  possible  complexities  in  such 
systems  came  from  the  work  of  Wang  and  Gomer[2]  who  found  remarkably  large 
effects  of  coadsorbed  oxygen  on  the  work  functions  and  desorption  temperatures 
of  Xe  on  W(110)  and  W ( 1 00 ) .  The  present  work  reveals  some  additional  surprises. 
In  brief,  it  was  found  that  desorption  could  be  followed  layer  by  layer  for  at 
least  3  ad-layers  and  that  in  each  case  desorption  was  initially  of  zero  order, 
until  roughly  0.75  of  a  given  layer  had  been  desorbed;  at  this  coverage  desorp¬ 
tion  became  first  order.  For  the  first  adsorbed  layer  a  very  small  amount  of 
preadsorbed  oxygen  sufficed  to  make  desorption  first  order  over  the  entire  Xe 
coverage  range.  Perhaps  even  more  surprising,  it  was  found  that  in  each  case 
the  activation  energy  of  desorption  was  substantially  less  than  the  correspond- 
i ng  enthalpy. 

EXPERIMENTAL 

The  apparatus  used  in  these  experiments  has  been  described  in  detail  pre¬ 
viously  [3].  Use  of  the  cryoshielded  effusion  source  made  it  possible  to  deposit 
Xe  only  on  the  front  surface  of  the  crystal,  which  was  kept  at  ~27K  during 
deposition.  Most  of  the  experiments  were  performed  with  the  cryoshield  slit  so 
wide  that  the  entire  length  of  the  crystal  received  a  gas  deposit.  fn  order  to 
avoid  possible  effects  from  the  ends  of  the  current  or  potential  leads  and  to 
improve  temperature  uniformity,  some  experiments  were  also  carried  out  with  a 
reduced  slit  length,  so  that  only  the  central  10  mm  portion  of  the  crystal, 
located  between  the  potential  leads,  received  a  Xe  deposit.  It  was  possible 
to  verify  this  adsorbate  geometry,  and  in  fact,  to  ascertain  the  correct  posi¬ 
tioning  of  the  crystal  in  front  of  the  source  by  means  of  Auger  measurements  in 
which  the  Xe  signal  was  measured  as  function  of  the  lateral  displacement  of  the 
crystal  relative  to  the  Auger  beam.  Temperature  programmed  and  isothermal 
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desorption  sequences  utilized  a  new  temperature  controller,  recently  described 
in  some  detail  [4].  For  most  experiments  desorption  rates  were  measured  directly 
with  a  UTI  C-100  quadrupole  mass  spectrometer.  For  coverages  of  <  I  monolayer 
the  amount  of  Xe  remaining  on  the  surface  was  also  determined  from  the  intensity 
of  the  Xe  MNN  Auger  transition.  In  these  measurements  heating  sequences  were 
interrupted  to  make  Auger  measurements  to  avoid  the  effects  of  magnetic  fields 
from  the  dc  heating  current  and  ESO  from  the  Auger  beam,  which  occurred  only  at 
T>70K. 

Relative  coverages  were  most  easily  determined  from  the  thermal  desorption 
spectra  themselves,  but  could  also  be  found  from  the  dosing  times,  once  that  for 
monolayer  coverage  had  been  established.  Oxygen  coverages  were  determined  from 
Auger  intensities,  relative  to  the  G/W  «  0.5  layer. 

The  temperature  calibration  of  the  crystal  has  been  described  previous  1 y [4] . 
RESULTS 

Fig.  I  shows  temperature  programmed  spectra  for  coverages  9,  2  <  6  ^  3* 
for  the  whole  front  surface  of  the  crystal.  3  distinct  peaks  are  clearly  dis¬ 
cernible,  with  those  corresponding  to  1  and  2  layers  remaining  unchanged,  while 
that  corresponding  to  the  3rd  layer  diminishes  in  size  with  decreasing  9^.  More¬ 
over,  the  leading  edge  Of  the  third  layer  peak  remains  unchanged,  indicating 
zero  order  desorption.  This  will  be  shown  in  greater  detail  by  the  isothermal 
results.  Fig.  2  shows  desorption  spectra  for  1  <9^2.  Here  the  spectrum  of 
the  first  layer  remains  unchanged,  while  that  of  the  second  decreases  in  size 
as  decreases.  The  leading  edge  of  the  second  layer  peak  again  remains  un¬ 
changed  as  02  decreases.  In  Fig.  1  the  time  scale  is  drawn  linearly  so  that  the 
areas  under  the  second  plus  third  layer  peaks  (for  a  full  3rd  layer)  are  twice 
that  of  the  first  layer.  In  Fig.  2  the  curves  have  been  redrawn  to  linear  tem¬ 
perature  scale  so  that  the  area  under  the  first  layer  peak  appears  too  small. 
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Fig.  3  shows  desorption  spectra  for  9  <  1.  Again  the  leading  edge  of  the  peak 
is  unchanged,  except  for  very  low  initial  coverages.  Essentially  identical 
results  were  obtained  by  dosing  only  the  central  portion  of  the  crystal. 

The  details  of  the  desorption  kinetics  are  brought  out  more  clearly  by 
isothermal  sequences.  Figs.  4-6  show  such  data  for  desorption  in  the  first, 
second  and  third  layers,  respectively.  In  each  case  the  desorption  rate  is 
almost  coverage  independent  for  9  >  0.7  (within  a  given  layer).  There  is  a 
sharp  break  at  this  critical  9,  and  for  lower  coverages  the  rate  -n  decreases 
rapidly  with  decreasing  9.  Figs.  7*9  show  that  plots  of  i^  vs.  time  (where 
iuc,  the  mass  spectrometer  current,  is  proportional  to  desorption  rate)  are 
then  linear,  indicating  first  order  desorption  in  this  regime.  These  results 
were  obtained  for  both  Xe  deposit  geometries.  Evidence  for  the  existence  of 
at  least  3  uniform  ad-layers  also  conges  from  the  fact  that  the  W  Auger  signal 
intensity  decreases  monoton i ca 1 1 y  with  Xe  coverage.  If  there  were  some  Xe 
crystallite  formation,  instead  of  adsorption  layer  by  layer,  this  would  not  be 
the  case. 

Even  small  amounts  of  preadsorbed  oxygen  lead  to  a  disappearance  of  the 
quasi-zero  order  regime  in  the  first  layer.  Fig.  10  shows  isothermal  desorption 
for  0/W  =*  0.5  (adsorbed  at  27  K  and  heated  for  30  sec  to  90  K  before  recooling 
and  Xe  adsorption)  for  a  saturated  Xe  monolayer.  Desorption  is  now  first  order 
over  the  entire  coverage  range.  Disappearance  of  the  quasi -zero  order  regime 
occurs  when  0/W  >  0.04. 

The  effect  of  small  amounts  of  oxygen  on  the  rate  suggests  that  the  zero 
order  regime  is  not  the  result  of  an  experimental  artifact,  connected  for  in¬ 
stance  with  pumping  speed  effects  or  with  temperature  anisotropies  over  the 
crystal.  The  first  of  these  points  could  be  checked  by  suddenly  turning  off 
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the  temperature  controller.  The  mass  spectrometer  signal  returned  to  zero  al¬ 
most  instantaneously.  The  second  point  was  checked  by  dosing  the  central  por¬ 
tion  of  the  crystal  to  monolayer  coverage  and  then  heating  to  66  K  for  90  sec, 
so  that,  according  to  curves  like  those  of  Fig.  A  the  remaining  coverage  should 
have  been  0.8.  The  crystal  was  then  displaced  parallel  to  its  major  dimension 
and  the  intensity  of  the  Auger  signal  determined  at  each  position.  No  difference 
could  be  discerned  within  experimental  error,  indicating  that  the  Xe  distribu¬ 
tion  remained  uniform.  It  is  most  unlikely  that  diffusion  could  be  rapid  enough 
at  60-80  K  to  achieve  a  substantial  redistribution  of  adsorbate.  Consequently, 
the  temperature  must  have  been  quite  uniform  over  the  central  10  mm  portion  of 
the  crystal . 

Experiments  were  also  carried  out  for  0  <  1  in  which  the  coverage  remain¬ 
ing  on  the  surface  was  determined  from  the  peak  to  peak  separation  of  the  second 
derivative  of  the  Auger  MNN  transition.  Results  obtained  in  this  way  are  shown 
in  Fig.  11  and  indicate  that  the  coverage  decreases  linearly  with  time  (as  it 
should)  in  the  quasi-zero  order  desorption  regime. 

In  order  co  obtain  the  zero-order  rate  constants  k  we  need  to  know  (c/A) 

o 

in  the  relation 


c  i 


MS 


A  5  k  A 
dt  o 


(1) 


where  is  the  mass  spectrometer  signal,  c  a  proport ional i ty  constant,  A  the 
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area  of  the  crystal  and  0  the  coverage  in  atoms/cm  .  It  is  known  from  the 
measurements  of  Wang  and  Gomer  [2]  that  the  maximum  monolayer  coverage  0O  * 

6.6  x  10^  at/cm^.  Since  the  total  area  under  curves  like  those  of  Fig.  3  cor¬ 


responds  to  Bo  we  can  find  (c/A)  from 

c/A  -  8  (t)/CiMS,t) 

where  0  ( t )  is  a  coverage  decrement,  obtained  from  the  relevant  fraction  of  the 


(2) 


total  area  under  the  i 


vs.  t  curve  corresponding  to  a  time  t.  Alternatively 


L* 


MS 

it  is  possible  to  determine  the  first  order  rate  constant,  kj  from  the  relation* 

(c/A)  T^s  *  -d9‘/dt'  3  k19Q  e  kl t  (3) 

i  i 

where  9q  is  the  coverage  where  first  order  desorption  starts  and  t  the  time 
measured  from  this  point: 


kj  ■  -  d ( 1 n  i MS) /dt ‘  (4) 

kj  can  then  be  used  in  Eq.  3, for  instance  at  t'  a  0  to  determine  (c/A).  Both 

methods  give  the  same  result  from  which  the  absolute  values  of  k  can  be  found. 

o 

Figs.  12-14  shows  plots  of  In  kQ  vs.  1 /T  for  the  initial  zero  order  desorption 
rates  in  the  first,  second,  and  third  layers.  Straight  lines  are  obtained  and 
yield  values  of  the  corresponding  activation  energies  which  are  listed  in  Table 
l .  Writing 


k  =  v  e 
o  o 


■E/kT 


(5) 


we  can  then  combine  these  activation  energies  with  the  absolute  values  of  k 

o 

to  obtain  the  correspond i ng  values,  also  listed  in  Table  1.  The  same  values 
were  obtained  by  plotting  the  logarithms  of  the  times  required  for  desorption 
to  some  constant  converage  (still  in  the  zero  order  regime)  vs.  l/T.  Figs. 
15*17  shows  plots  of  In  k^  vs.  l/T  for  the  first  order  desorption  regimes  in 
the  first,  second,  and  third  layers  respectively.  The  corresponding  activation 
energies  and  frequency  factors  are  listed  in  Table  1. 


As  already  pointed  out,  preadsorbed  oxygen  in  amounts  G/W  >  0.04  led  to 
first  order  desorption  in  the  first  layer.  Fig.  18  shows  a  temperature 
programmed  desorption  spectrum,  starting  with  2.7  monolayers  of  Xe  for  sur¬ 
faces  precovered  with  0/W  =  0.5  preheated  to  90  K  and  900  K.  For  both  con¬ 
ditions,  the  desorption  peak  of  the  first  layer  is  shifted  to  increased 
temperatures.  For  0/W  =  0.5*  900  K,  there  is  also  a  long,  low  temperature 
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tail,  which  will  be  discussed  presently.  It  is  interesting  that  adsorbed  oxygen 
also  affects  the  second  Xe  layer,  but  not  the  third.  Fig.  18  shows  clearly 
that  the  second  layer  peak  is  shifted  to  lower  temperature  by  the  presence  of 
oxygen,  while  the  third  layer  peak  is  not  affected.  Oxygen  precoverage  does 
not  remove  the  zero  order  regime  for  the  second  layer. 

Activation  energies  and  prefactors  for  various  oxygen  precoverages  heated 
to  90  K  were  obtained  for  the  first  Xe  layer.  The  resu 1 ts, (Fi gs .  19  and  20), 
show  considerable  scatter, but  indicate  a  slight  increase  in  E  and  a  decrease  in 
v  with  increasing  coverage.  Perhaps  more  significant  are  the  values  of  k^  at 
a  fixed  temperature  of  76.5  K  vs.  oxygen  coverage,  shown  in  Fig.  21.  These 
indicate  a  systematic  decrease  of  in  the  range  0.05  <0/W<0.12,  followed  by 
a  1 evel 1 i ng  off. 

Isothermal  desorption  from  surfaces  covered  with  oxygen  to  the  extent  of 
0/W  =  0.5  preheated  to  900  K  reveal  three  distinct  desorption  regimes.  Starting 
with  a  full  Xe  monolayer  on  the  central  portion  of  the  crystal , there  is  desorp¬ 
tion  of  0.25  layers  in  the  temperature  range  62  -  75  K.  The  desorption  of  this 

fraction  can  be  approximated  by  2  first  order  regimes  with  activation  energies 
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of  2.8  and  2.7  Kcal/mol  and  pre-exponentials  of  5  x  10'  and  2  x  10  sec  re¬ 
spectively.  When  the  entire  face  of  the  crystal  was  dosed, the  results  were 
E  *  3-5  and  3.4  Kcal/mole  respectively  and  v  =  2  x  10^°  for  both.  After  desorp¬ 
tion  of  this  f ract ion  ,  (for  central  dos i ng) , another  0.2  layers  desorb  in  the 
range  72-78  K,  with  E  =  5.8  Kcal/mole  and  v  =  6  x  10^5  sec  \  Finally,  the  re¬ 
mainder,  corresponding  to  0.55  layers , desorbs  in  the  range  78-86  K  with  E  =  6.15 
Kcal/mole  and  v  »  2  x  10^  sec  ^  .  When  the  entire  front  surface  of  the  crystal 
was  dosed,  these  last  2  states  were  not  resolved  and  a  single  regime  with  E  = 

6.4  Kcal/mole  and  v  *  8  x  lO^sec  ^  was  seen  for  the  terminal  0.55  monolayers. 
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This  complex  behavior  is  indicated  qualitatively  by  the  relevant  desorption  spec¬ 
trum  of  Fig.  18.  Fig.  22  shows  plots  of  In  kj  vs.  1 /T  for  the  three  regimes, 
including  the  decomposition  of  the  lowest  temperature  regime  into  2  exponentials. 


The  results  for  precoverage  with  oxygen  at  90  and  900  K  are  in  remarkably 
good  agreement  with  the  qualitative  results  obtained  by  Wang  and  Gomer  [2]  for 
similar  conditions,  although  these  authors  could  only  obta i n. somewhat  indirectly, 
crude  “step  desorption11  spectra.  In  particular,  they  concluded  that  heating  an 
oxygen  layer  to  900  K  led  to  3  distinct  adsorption  states  for  Xe,  and  that  pre¬ 
heating  to  90  K  increased  the  desorption  temperature  of  Xe. 

It  is  interesting  that  the  desorption  data  obtained  for  the  first  order 
regime,with  only  the  center  of  the  crystal  covered  with  Xe, yield  more  finely 
resolved  regimes  probably  because  of  the  slightly  more  uniform  temperature  near 
the  center  of  the  crystal  and  the  absence  of  surface  heterogeneities  like  the 
ends  of  leads. 


DISCUSSION 


Zero-Order  Desorption 


The  most  striking  feature  of  the  present  results  is  the  occurrence  of  zero- 
order  desorption,  followed  by  a  sharp  transition  to  first  order  desorption. 
Zero-order  desorption  has  also  been  seen  in  a  number  of  other  systems,  e.g., 

Hg  on  W(100)  [5],  Ag  on  Ge  (111)  [6],  Zn  on  GaAs  [7],  and  Xe  on  graphite  (1000) 
[8].  Several  explanations  to  account  for  this  behavior  have  been  advanced,  hut 
none  of  these  really  explain  our  results  satisfactorily. 

We  will  presently  show  that  the  formation  of  close  packed  islands  is  prob¬ 
ably  involved,  but  not  in  a  simple  way.  We  first  indicate  why  most  island  mecha¬ 


nisms  cannot  explain  the  observed  behavior. 

1)  Desorption  from  the  edge  of  islands,  the  number  of  islands  N  remaining 
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fixed.  In  this  case  desorption  is  of  1/2  order  and  can  be  shown  to  have  the 
form 

dg/dt  »  -  =»  2  kj  (tt/9j  )  ^Nn^2,  (6) 

where  n  is  the  total  number  of  atoms  on  the  surface,  and  g  the  number  of  de- 

2 

sorbed  atoms.  9^  is  the  (constant)  density  in  atoms/cm  within  each  island, 
k^  is  the  desorption  rate  constant. 

2)  Desorption  from  the  second  1 ayer 9 ( re  1  at i ve  to  that  from  which  desorp¬ 
tion  occurs);  i.e.  from  the  tops  of  islands,  N  =  constant.  If  the  second  layer 
is  populated  from  the  edges  of  islands,  the  kinetics  are  as  in  (1).  If  the 
second  layer  is  populated  from  within  the  island,  kinetics  are  first  order,  as 
can  be  shown  trivially. 

3)  Detachment  into  a  dilute  phase,  desorption  from  the  dilute  phase, 

N  =  constant. 

k]  k2 

Dense  Islands  -*■  Dilute  Phase  Gas  Phase 

a)  ;  (2) 

-i 


We  have  for  this  assumed  mechanism 

dn2/dt  **  kj  N 

-dn^/dt  -  k_j  N2ur  (n2/A2)  +  k^n2 

where  r  is  the  mean  island  radius  and  A2  the  area  of  the  dilute  phase. 

steady  state  for  the  total  number  of  atoms  in  the  dilute  phase,  n2 

k.  N27ir  A„ 

n~  =  _\ _ 2 _ 

kH  N2ttt  +  k2  A2 


(7) 

(8) 

Assuming 

(9) 


Since 


r  =  (A  1/2 

'e.TT 


do) 


where  is  the  number  of  atoms  in  an  average  island,  and  9  the  (constant) 
density  within  each  island,  we  obtain 

da.  dn  k.  ki  «<-/9,i'/2",l/2  (*-»,) 

dt  dt  Id  ml  /«  ’I  1/2  1/2 

k-l  M2('/01)  "l  •+k2(A-A1) 


(ID 
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where  we  have  replaced  by  A-Aj ,  A^  being  the  area  occupied  by  islands.  If 

k_]  N2C/9,)'72  n^2  «k2  (A-A1 )  (12) 

we  obtain  Eq.  (6),  although  in  principle  now  refers  to  detachment  only. 


If  inequality  (12)  is  reversed  we  find 

-  ST"  K  k2  (A-V 


(I!) 


where 


^1-^-9 

k  A  02 

-1  t 


OM 


with  the  second  equality  following  from  Eq.  (9)  with  the  reversed  inequality 
(12).  Under  these  conditions  9^  -  constant.  Since 

-dn  =  -6^dA^  -92  d  A2  05) 

and  since 

dA}  =  -  dA2  (16) 

we  easily  find  from  Eq.  (13)  that 

Kk 

‘  HT  =  S  =  K  k2  V°)  exp  (17e2t}  (17) 

i.e.  exponentially  increasing  desorption,  assuming  a  small  but  finite  area 
A^  (0)  available  to  the  dilute  phase  at  time  t  =  0.  When  only  the  dilute  phase 
remains  desorption  would  become  first  order,  i.e.  dg/dt  would  then  decrease 
exponentially.  Thus,  none  of  the  cases  so  far  considered  lead  to  anything 
resembling  zero  order  kinetics.  While  it  is  true  that  Eq.  13  corresponds  to 
constant,  equilibrium  concentration  in  t m  dilute  phase,  the  area  of  the  dilute 
phase  increases  as  desorption  proceeds  thus  leading  to  exponentially  increasing 
desorption  as  just  shown. 

4)  Attractive  interactions  between  ad-atoms  without  explicit  island 


f 


i 


format i on . 
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In  this  case  the  qualitative  behavior  is  in  the  right  direction.  As  the 
adsorbate  concentration  decreases,  the  ad-ad  interactions  play  a  decreasing 
role  thereby  decreasing  ^cjes  and  compensating  for  the  decrease  in  rate  which 
would  otherwise  occur.  Figs.  2 3 “24  show  the  results  of  Monte-Carlo  calculations 
for  various  assumed  Xe-Xe  interaction  parameters.  If  the  interaction  is  too  weak, 
quasi-first  order  desorption  results.  If  it  is  too  strong, an  initial  increase 
in  rate  followed  by  a  decrease  occurs.  Such  peaking  was  also  found  by  Jones 
and  Perry[9j.  In  no  instance  was  it  possible  for  us  to  obtain  behavior  in  which 
the  change  from  quasi-zero  order  to  first  order  was  either  abrupt  or  occurred 
at  the  correct  coverage.  The  results  were  not  improved  by  including  diffusion 
either  with  the  same  or  with  different  interaction  parameters  as  for  desorption. 
Some  details  of  these  calculations  are  given  in  Appendix  I. 

We  come  finally  to  two  mechanisms  involving  islands  which  can  explain  the 
results.  Desorption  from  the  edges  of  islands,  either  directly  into  the  gas 
phase^r  first  into  a  dilute  phase  with  high  probability  of  subsequent  desorp¬ 
tion,  leads  to  zero  order  kinetics  if  the  total  periphery  of  all  islands  remains 
constant.  This  is, in  fact, not  unreasonable.  At  a  critical  coverage,  depending 
slightly  on  temperature,  islands  break  up  completely  and  desorption  then  proceeds 
via  first  order  kinetics.  This  transition  can  be  quite  sharp  accounting  for  the 
sharp  break  in  the  curves  of  Figs.  A-6.  In  the  zero  order  regime  we  have  for 
th i s  mechan i sm 

dg/dt  =  kQ  A  =  2tt  rNkj  =  27rr  (N/A)Ak]  (18) 

where  all  symbols  have  the  same  meaning  as  in  Equations  6-17*  Thus, 

k  =  2 nr  (N/A)k. 
o  1 

As  here  defined  k^  is  the  rate  of  evaporation  per  unit  line  length.  In  terms 
of  the  effective  frequency  of  atomic  vibration,  Vj ,  and  the  activation  energy, 
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E  per  a tom, we  can  write 

k,  =  (v,  e'E/kT)/dXe  (19) 

-8 

where  A  3  x  10  cm  represents  the  perimeter  length  per  Xe  atom.  Thus,  we 
can  write 

Vo  =  ^TTr  ^l^Xe^  (20) 

We  assume,  to  keep  the  kinetics  zero  order,  that  2tt r  (N/A)  *  constant. 

12  “1 

If  we  assume  that  -  10  sec  we  see  from  Eq.  20  and  from  the  experimental 

25  r  _  i 

value  \>  *  1  x  10  for  the  first  layer  that  2:rr  (N/A)  =  3  x  10pcm  .  Since 

o 

(N/A)7rr2  (21) 

we  find  that  r  =  6  x  10  ^cm.  Thus,  the  assumption  of  desorption  from  island 
max 

edges  with  the  total  periphery  constant  is  consistent  with  the  experimental 
value  of  and  a  reasonable  choice  for  if  the  maximum  average  island  radius 
is  "JO  ^cm.  For  the  second  and  third  layers, even  larger  values  of  r  are 
a  1 1  owed . 

There  is  a  second  possible  mechanism  involving  islands.  It  is  based  on  the 
following  assumptions: 

1)  Desorption  occurs  only  when  a  Xe  atom  has  become  detached  from  an  island, 
i.e.,  is  in  the  dilute  phase.  2)  At  the  desorption  temperature ,Xe  atoms  move 
so  rapidly,  and  their  concentration  in  the  dilute  phase  is  so  small  that  their 
mean  free  path  is  limited  by  collisions  with  islands  to  which  they  become 
reattached  upon  collision.  The  decrease,  if  any,  in  detachment  rate  because  c " 
decreasing  island  perimeter  is  then  exactly  balanced  by  an  increase  in  the  mean 
time  between  atom-island  collisions  since  the  former  is  proportional  to  island 
radius  and  the  latter  inversely  proportional  to  this  quantity.  The  time  t 
between  collisions,  i.e.,  the  mean  free  time  of  Xe  atom  in  the  dilute  phase  is: 
tc  -  [  V  (d  +  dXe)  (N/A)] 


(22) 
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where  d  is  the  mean  island  diameter  and  dv  =  3  x  10  cm  that  of  a  Xe  atom. 

Xe 

v  is  the  mean  tnermal  velocity  of  Xe  atoms  in  the  dilute  phase.  The  desorp¬ 
tion  rate  is  then 

,  t  71  k,  k.  A 

=  —  T  d-Nk  «  — - - -  (  -  -  ^ - )  (21) 

dt  1  -  K  d  +  dv  '  Uj5; 

t  v  Xe 

where  t  i  l/k^  is  the  lifetime  of  a  Xe  atom  in  the  dilute  phase  and  kj  the  de¬ 
tachment  rate  constant  from  the  island  edge  into  the  dilute  phase.  If  d^e<<  d 
desorption  is  zero  order  without  any  assumptions  about  the  constancy  of  the 
island  perimeter  or  changes  in  the  number  of  islands.  A  sharp  transition  to 
first  order  will  occur  if  islands  disappear  while  d  is  still  somewhat  larger 

than  dv  .  For  this  mechanism 
Xe 

k0  =  kl  k2~/  7  (2A) 

and  the  experimental  activation  energy  is  the  sum  of  the  detachment  and  evapora¬ 
tion  activation  energies.  Consequently,  we  may  write 

vo  =  V1  v2^  ^  x  10  (25) 

_  1/2  3 

Assuming  v  =  (kT/2m)  =  A. 7  x  1  Cr  cm/sec  at  70  K  and  taking  from  the 

1 2  - 1  8-1 

first  order  regime  to  be  A  x  10  sec  ,  we  find  =  1  x  10  sec  for  this 
mechanism.  This  is  rather  small,  but  not  impossible. 

It  is  interesting  that  the  activation  energies  for  zero  and  first  order 
kinetics  are  almost  identical  within  experimental  error,  and  that  the  presence 
of  preadsorbed  0  also  affects  E  only  very  slightly.  This  indicates  that  very 
small  energy  differences  determine  which  regime  will  be  observed,  and  also  ex¬ 
plains  why  zero  order  desorption  can  be  observed, even  at  full  coverage* where 
presumably  the  island  perimeters, (or  the  distance  between  islands)*are  zero 
not  counting  the  external  edge  of  the  layer.  Adsorbed  oxygen  presumably  in¬ 
hibits  island  formation,  possibly  because  theO-Xe  interaction  predominates 
over  the  Xe-Xe  interaction.  Although  the  Monte-Carlo  calculations  mentioned 
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earlier  do  not  give  a  quantitative  description  of  the  observed  desorption  be¬ 
havior,  they  come  rather  close.  It  is  probably  significant  that  the  best  fit 
is  obtained  for  Xe-Xe  interactions  of  45  cal/mole,  a  very  small  energy  relative 
to  the  (4000  cal/mole)  activation  energy  of  desorption  in  the  first  layer. 

Activation  Energies  of  Desorption 

The  enthalpy  of  sublimation  of  bulk  Xe  is  3*7  Kcal/mole.  Thus,  even  in  the 
first  layer, the  activation  energy  is  only  very  slightly  greater  than  the  Xe  sub¬ 
limation  enthalpy.  In  the  second  and  third  layers, the  activation  energies  of  3*2 
and  2.9  Kcal  are  substantially  less.  Since  the  entropies  of  bulk  Xe  and  adsorbed 
Xe  layers  must  be  comparable,  it  follows  that  the  enthalpies  of  adsorption  must 
exceed  that  of  bulk  Xe  or  the  adlayers  would  be  unstable  with  respect  to  the 
latter.  There  is  absolutely  no  evidence  of  such  instability  up  to  at  least 
three  adsorbed  layers.  Consequently,  we  can  conclude  that  the  activation  ener¬ 
gies  in  the  second  and  third  layers  are  definitely  less  than  the  correspond i ng 
enthalpies.  Even  in  the  first  layer  this  is  almost  certainly  the  case.  Pre¬ 
vious  work[10] ,  suggests  that  the  enthalpy  of  adsorption  in  the  first  layer  is 
closer  to  10  Kcal/mole  than  to  4  Kcal/mole. 

Thus,  we  have  here  a  quite  clearcut  case  of  activation  energies  less  than 
corresponding  enthalpies.  A  similar  situation  was  encountered  by  us  for  weakly 
adsorbed  molecular  oxygen  both,  oi-O^  and  physisorbed  0^,  which  exist  on  a  W ( 1 10) 
surface  precovered  with  atomic  0[4] .  Such  behavior  has  been  predicted  theoreti¬ 
cally  by  Freed  and  coworkers  [11].  Their  explanation,  at  least  qualitatively, 
is  the  following.  The  potential  energy  curve  cor respond i ng  to  the  motion 
of  the  adsorbate  normal  to  the  surface  deviates  from  that  of  a  harmonic 
oscillator  as  the  surface-adsorbate  distance  increases.  In  the  enharmonic 
region  the  density  of  states  increases  rapidly  with  increasing  energy 
and  if  this  degree  of  freedom  is  excited  to  a  certain  critical  energy 


which  may  lie  quite  appreciably  below  that  corresponding  to  desorption,  the 
probability  that  the  adsorbate  will  gain  enough  energy  from  multi -phonon  pro¬ 
cesses  to  desorb  is  much  greater  than  the  probability  that  it  will  lose  energy 
and  return  to  the  ground  state  by  emitting  one  phonon  at  a  time.  Thus,  the 
activation  energy  measured  in  a  kinetic  experiment  corresponds  to  this  critical 
energy  rather  than  to  the  enthalpy  of  desorption.  The  extreme  case  of  the 
situation  just  described  can  be  represented  as  a  two  step  mechanism  for  desorp¬ 
tion:  step  1  takes  the  system  to  the  critical  energy,  and  step  2  corresponds  to 
desorption  from  this  point.  If  step  1  is  much  slower  than  step  2,the  overall 
rate  is  governed  by  its  rate. 

This  can  be  formalized  by  the  requirement  that 


-Ej/kT 


-E^/kT 


(26) 


or,  since  +  E^  =  H,  where  H  is  the  enthalpy  of  desorption,  that 

,  _  (2E,-H)/kT  ,97x 

Vj/v^  «  e  1  (27) 

Effect  of  Oxygen  on  Xe  Desorption 

Preadsorbed  oxygen  heated  to  90  K  appears  to  have  only  a  small  effect  on 
the  activation  energy  of  Xe  desorption.  in  view  of  the  preceding  discussion, 
it  is  not  possible  to  equate  this  with  an  equally  small  effect  on  the  enthalpy, 
but  this  seems  probable.  Discussions  of  why  adsorbed  0  could  affect  the  enthalpy 
of  Xe  adsorption  have  been  given  prev ious 1 y [2]  and  we  have  little  to  add. 

The  situation  for  0  layers  heated  to  900  K  is  more  intriguing.  Approximately 


0.45  of  the  Xe  layer  now  desorbs  with  a  smaller  activation  energy  than  on  clean 
W,  while  "0-55  desorbs  with  a  higher  one.  For  this  fraction^the  high  value  of 
v  suggests  that  the  activation  energy  is  now  fairly  close  to  the  enthalpy  in 


view  of  Eqs.  (26)  and  (27).  While  no  real  conclusions  about  the  effect  on  the 
enthalpy  can  be  drawn,  it  appears  that  heating  to  900  K  leads  to  surface  recon¬ 
struction  which  permits  a  part  of  the  adsorbate  to  be  in  such  contact  with  W 
atoms  as  to  increase  its  binding  energy  while  forcing  a  smaller  fraction  of  Xe 
to  be  adsorbed  in  a  way  which  actually  decreases  its  binding  energy.  It  is 
planned  to  investigate  these  points  by  means  of  UPS  and  XPS. 

CONCLUSION 

The  results  discussed  here  indicate  that  the  adsorption  of  Xe  on  W ( 1 10) 
represents  a  very  complicated  system.  The  principal  findings  are  the  occurrence 
of  zero  order  desorption  in  each  adsorbed  layer  and  clear  evidence  for  activa¬ 
tion  energies  less  than  the  cor respond i ng  enthalpies.  The  latter  finding  indi¬ 
cates  that  one  must  be  very  cautious  about  equating  activation  energies  with 
enthalpies  of  desorption,  particularly  for  weak  adsorption. 
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A  Monte-Carlo  simulation  of  desorption  was  carried  out  as  follows:  A  square 
lattice  was  assumed  with  periodic  boundary  conditions.  Pairwise  nearest  neigh¬ 
bor  interactions  J  were  assumed.  Full  coverage  at  the  start  of  each  run  was 
chosen.  In  each  complete  computer  pass  each  site  was  interrogated  in  random 
sequence  once  and  only  once.  The  desorption  probability  for  each  site  in  one 
such  interrogation^  (i .e. ,  per  computer  step)?was  chosen  to  be 

W  =  t  ^  exp  (nJ/kT)  for  an  occupied  site  (Al) 

des 

W  =  0  for  an  empty  site 

Attractive  nearest  neighbor  interactions  correspond  to  negative  values  of  J. 

corresponds  to  the  mean  lifetime  with  respect  to  desorption  in  the  absence 
of  interactions,  i.e.  at  zero  coverage: 

Tdes  *  v  exP  ^des^"^  (in  seconds)  (^) 

where  visa  frequency  factor  and  E^es  the  zero  coverage  activation  energy  of 
desorption.  Since  es  is  expressed  in  real  time,  this  fixes  the  scale  of  com¬ 
puter  time  as  1  second  per  step.  This  assignment  is  adequate  for  x  >>  1  second, 
as  was  the  case  here.  (in  any  case,  a  constant  factor  converting  r  from  real 
to  an  arbitrary  time  scale  could  have  been  used  if  this  condition  were  not  met.) 
The  probability  of  desorption  according  to  Eq.  Al  was  then  determined  in  a 
standard  manner  by  assigning  the  appropriate  weight  to  W  in  the  interval  0-1 
and  picking  a  random  number  in  this  interval  by  means  of  a  DEC  pseudo-random 
number  generator.  The  results  of  many  runs  from  full  coverage  to  zero  coverage 
for  a  20  x  20  lattice  were  averaged  to  obtain  the  data  shown  in  Fig.  23  whicn 
also  shows  the  results  of  a  real  experiment.  The  Monte-Carlo  results  for 
J  =  -  0.0**  Kcal/mole  come  closest  to  simulating  the  real  experiment,  but  lack 
the  sharp  transitions  from  zero  to  first  order. 


Additional  Monte-Carlo  experiments  were  carried  out  to  include  the  effect 
of  diffusion.  In  pr inci pie, thi s  can  be  done  by  sampling  each  site  as  before 
and  dividing  the  number  interval  0  -  1  into  3  segments  corresponding  to  the 
probability  of  diffusion,  desorpt ion,or  of  no  event.  The  probability  of  desorp¬ 
tion  is  then  given  by  Eq.  A1 ,  that  of  diffusion  is  given  by 

Wdiff  *  Tdiff  exp^n  ^dif^k^)  if  the  site  is  occupied 

=0  if  the  site  is  empty  (A3) 

If  the  random  number  generator  signals  diffusion,  a  nearest  neighbor  site  is 
then  picked  at  random  and  if  this  site  is  empty  ,d i ff us  ion  to  it  is  carried 
out.  If  the  site  to  which  diffusion  is  proposed  is  filled,  no  event  occurs. 

The  difficulty  with  this  approach  is  that  the  value  of  is  so  much  less 

than  that  of  that  each  unit  step  can  no  longer  be  chosen  as  1  second.  If 

the  values  of  Chen  and  Gomer[12]  for  Xe  diffusion  are  chosen  =  10^  sec  ^ 

-ft  2  - 1  L 

(assuming  D  =  10‘°cm  sec  0)  and  E  , .  r  «  1  Kcal/mole,  =  1  x  lO^sec. 

o  a  1 1  □  1 1 

2  - 1 

Experiments  with  an  unrealistically  small  value  of  =  6.25x10  sec  and 

- 1  3  2  - 1 

^dif  *  *  Kcal/mole  (giving  0q  =  5  x  10  J  cm  sec  ,  and  *  2)  were  in  fact 

carried  out,  and  led  to  very  little  change  from  the  previous  results.  Note  that 
here  =  0- 1  »  so  that  quasi  -equi  1  i  brat  ion  occurs  between  desorption  events. 

A  more  realistic  approach,  given  the  disparity  between  and  con¬ 

sists  of  making  one  desorption  step  as  before,  allowing  the  sample  to  equili¬ 
brate  by  means  of  standard  Monte-Carlo  procedure ,  (not  counting  in  any  way  the 
time  for  this  equ i 1 i brat  ion),  and  then  proceeding  to  the  next  desorption  step, 
etc.  The  results  obtained  in  this  way  were  not  significantly  different  from 
the  previous  ones. 

We  would  like  to  acknowledge  many  helpful  discussions  with  Stanley  Brooks 


and  Jayanth  Banaver  about  the  Monte  Carlo  calculations. 
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FIGURE  CAPTIONS 


1)  Temperature  programmed  desorption  for  initial  coverages  in  excess  of  2  mono- 
layers.  Note  the  constant  leading  edge  of  the  desorption  peak  of  the  third 
layer.  The  variation  in  peak  height  for  the  second  layer  peak  is  probably 
due  to  some  overlap  from  the  third  layer  peak.  Entire  front  face  of  crystal 
covered  with  Xe. 

2)  Temperature  programmed  desorption  for  coverages  in  excess  of  1  layer  but  less 
than  2  layers.  Note  the  constant  leading  edge  of  the  second  layer  peak.  The 
apparent  difference  in  area  under  the  first  and  full  second  layer  peaks  is 
discussed  in  the  text.  Entire  front  face  of  crystal  covered  with  Xe. 

3)  Temperature  programmed  desorption  for  9  ,<  1  monolayer.  Note  the  constant 
leading  edge  of  the  peaks.  Entire  front  face  of  crystal  covered  with  Xe. 

A)  Isothermal  desorption  for  initial  coverage  of  1  monolayer.  Note  the  almost 
constant  rate  of  desorption  to  9  3  0.3*  Initial  fast  decay  is  the  result  of 
temperature  overshoot.  Entire  front  face  of  crystal  covered  with  Xe. 

5)  Isothermal  desorption  within  the  second  layer.  Entire  front  face  of  crystal 
covered  with  Xe. 

6)  Isothermal  desorption  within  the  third  layer.  Entire  front  face  of  crystal 
covered  with  Xe. 

7)  Sem* logar i thmi c  plot  of  mass  spectrometer  signal  vs.  time  for  first  order 
regime  for  first  Xe  layer.  Only  the  central  portion  of  front  face  of  crystal 
covered  with  Xe. 

8)  Semi logr i thmi c  plot  of  mass  spectrometer  signal  vs.  time  for  first  order 
regime  in  the  second  layer.  Entire  front  face  of  crystal  covered  with  Xe. 

9)  Semi logr i thmi c  plot  of  mass  spectrometer  signal  vs.  time  for  first  order 
regime  in  the  third  layer.  Entire  front  face  of  crystal  covered  with  Xe . 

10)  Mass  spectrometer  signal  vs.  time  for  isothermal  desorption  of  Xe  in  the 
first  layer  for  a  surface  precovered  with  oxygen  and  heated  to  90  K  before 
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Xe  adsorption.  Entire  front  face  of  crystal  covered  with  oxygen  and  Xe. 

Ratio  of  Xe  to  W  Auger  intensities,  determined  from  peak  to  peak  amplitudes 
for  the  zero  order  regime  in  the  first  layer. 

Plot  of  In  for  the  first  layer  vs.  l/T.  All  open  points  refer  to  runs  with 

the  entire  front  face  of  the  crystal  dosed,  full  points  to  a  Xe  deposit  only 
on  the  central  part  of  this  face. 

Plot  of  In  kQ  for  the  second  layer  vs.  l/T.  All  open  points  refer  to  runs 
with  the  entire  front  face  of  the  crystal  dosed,  full  points  to  a  Xe  deposit 
only  on  the  central  part  of  this  face. 

Plot  of  In  for  the  third  layer  vs.  l/T.  All  points  were  taken  with  the 

entire  front  face  of  the  crystal  dosed  with  Xe. 

Plot  of  In  kj  vs.  l/T  for  the  first  layer.  Open  points  and  solid  line  refer 
to  entire  front  face  of  crystal  covered  with  Xe,  solid  points  and  dashed  line 
to  Xe  deposited  only  on  the  central  portion  of  the  face. 

Plot  of  In  kj  vs.  l/T  for  the  second  layer.  Open  point  and  solid  line  refer 
to  entire  front  face  of  crystal  covered  with  Xe ,  solid  points  and  dashed  line 
to  Xe  deposited  only  on  the  central  portion  of  the  face. 

Plot  of  In  kj  vs.  l/T  for  the  third  layer.  Entire  front  face  of  crystal  dosed 
with  Xe. 

Temperature  programmed  desorption  for  2.7  monolayers  of  Xe  on  surfaces  pre¬ 
covered  with  oxygen  and  heated  as  indicated  on  the  figure.  Oxygen  coverages 
in  all  cases  0/W  3  0.5.  Note  the  effects  on  both  the  first  and  second  layer 
peaks . 

Activation  energy  for  first  order  regime  in  the  first  layer  as  function  of 
oxygen  coverage,  expressed  as  0/W,  for  the  oxygen  layer  heated  to  90  K  before 
Xe  adsorption.  Entire  front  face  of  crystal  covered  with  oxygen  and  Xe. 


FIGURE  CAPTIONS 


20)  Prefactor  of  first  order  rate  constant  in  the  first  layer  as  a  function  of 

oxygen  coverage  expressed  as  0/W,  for  the  oxygen  layer  heated  to  90  K  before 

Xe  adsorption.  Ent i re  front  face  of  crystal  covered  with  oxygen  and  Xe. 

21)  In  at  76.5  for  first  order  desorption  in  the  first  layer  as  a  function  of 

oxygen  coverage  expressed  as  0/W,  for  the  oxygen  layer  heated  to  90  K  before 

Xe  adsorption.  Entire  front  face  of  crystal  covered  with  oxygen  and  Xe. 

22)  In  vs.  1/T  for  desorption  in  the  first  layer  from  a  surface  precovered 
with  oxygen,  0/W  =  0.5  and  heated  to  900  K  before  Xe  adsorption.  Note  3 
regimes.  Entire  front  face  of  crystal  covered  with  oxygen  and  Xe. 

Open  points  and  solid  lines  refer  to  entire  front  face  of  crystal  covered 
with  Xe,  solid  points  and  dashed  lines  to  Xe  deposited  only  on  the  central 
portion  of  the  face. 

23)  Monte-Carlo  calculations  of  desorption  rate  An/At  in  atoms/second  vs.  real 

time  for  various  assumed  interaction  energies  J,  shown  on  figure.  The  zero 
coverage  desorption  mean  lifetime  *  23.1  seconds.  An  experimental 

curve  is  also  shown. 

2k)  Monte-Carlo  calculations  of  desorption  for  various  assumed  adsorbate-adsorb¬ 
ate  interactions  in  desorption  and  different  interactions  in  diffusion  J, 

d 

shown  on  figure,  x,  =  23.1  sec,  x..r  «  2.1  sec. 
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